Bernoulli 14(4), 2008, 1134-1149 
DOI: 10.3150/08-BEJ139 



The Dagum family of isotropic correlation 
functions 

CHRISTIAN BERG 1 , JORGE MATEU 2 '* and EMILIO PORCU 1 ' 2 " 

1 Department of Mathematics, University of Copenhagen, Universitetsparken 5, DK-2100 Copen- 
hagen, Denmark. E-mail: Berg@math.ku.dk 

2 Department of Mathematics, Universitat Jaume I, Campus Riu Sec, E-12071 Castellon, Spain. 
E-mail: mateu@mat.uji.es; porcu@mat.uji.es 

A function p : [0, oo) — > (0, 1] is a completely monotonic function if and only if p(||x|| 2 ) is positive 
definite on R d for all d and thus it represents the correlation function of a weakly stationary 
and isotropic Gaussian random field. Radial positive definite functions are also of importance 
as they represent characteristic functions of spherically symmetric probability distributions. In 
this paper, we analyze the function 



p((3,y)(x) = l- x>0, /3, 7 >0, 



called the Dagum function, and show those ranges for which this function is completely mono- 
tonic, that is, positive definite, on any d-dimensional Euclidean space. Important relations arise 
with other families of completely monotonic and logarithmically completely monotonic func- 
tions. 

Keywords: Bernstein function; completely monotonic function; Dagum family; isotropy; 
logarithmically completely monotonic function; Stieltjes transform 

1. Introduction 

Positive definite functions on groups or semigroups have a long history and arc present 
in many applications in probability theory, operator theory, potential theory, moments 
problems and several other areas. They enter as an important chapter in all treatments 
of harmonic analysis, and can be traced back to papers by Caratheodory, Herglotz, 
Bernstein and Matthias, culminating in Bochner's theorem from 1932-1933. See [24] for 
details. 

Schoenberg's theorem ([25]) explains the possibility of constructing rotation- invariant 
positive definite and conditionally negative definite functions on Euclidean spaces via 
completely monotonic functions and Bernstein functions. Positive definite radial functions 
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play a very important role in probability theory and statistics. They occur both as 
the characteristic functions of spherically symmetric probability distributions ([8], and 
references therein) and as the correlation functions of isotropic random fields on Euclidean 
spaces ([18]). 

Among the valuable contributions to the development of random field theory at- 
tributable to eminent scientists, a crucial part is recovered by those regarding the con- 
struction of covariance/correlation functions associated with (wide sense) stationary ran- 
dom fields (RF for short). Covariance functions are fundamental in describing the depen- 
dence between the values of a spatial RF at locations separated by a lag vector x e M. d . 
The (very often) complex spatial interactions and types of dependencies call for sim- 
plifying assumptions and, typically, the spatial phenomenon is considered a realization 
of a stationary and additionally isotropic RF so that correlations between observations 
separated by some vector x G K d become a function of only the Euclidean distance ||x|| . 

Among the correlation models proposed in the literature, two are of particular impor- 
tance for RF theory: the Matern class ([16]) and the generalized Cauchy class ([12], and 
references therein). 

[27] argues convincingly that the Matern covariance function possesses desirable mathe- 
matical features that make this function very appealing to describe spatial dependencies. 
As for the generalized Cauchy class, it has been recently revisited by [13], who show 
heuristically that this class allows for the decoupling of the local and global behavior of 
the associated weakly stationary and isotropic RF. The Matern and Cauchy classes are 
permissible for any dimension of the associated Euclidean space M. d , d>l, and only a 
few other models are available in the literature: among them, it is worth citing the Gaus- 
sian model and the exponential model, associated with the so-called Ornstein-Uhlenbeck 
process or continuous autoregressive of the first order, that are additionally special cases 
of the Matern class for particular settings of the parameters indexing this correlation 
function. 

Recently, [23] introduced the Dagum class of radially symmetric correlation functions. 
Roughly speaking, for (3 > 0, 7 > 0, they considered the function 



and showed permissibility criteria for this function in a three-dimensional Euclidean 
space. Next, [17] showed that the Dagum class allows the fractal dimension and the Hurst 
effect of the associated weakly stationary Gaussian RF to be treated independently, by 
using the procedure suggested in [13]. Comparisons between Dagum and Cauchy models 
are considered in Section 5. 

It is well known that correlation functions must be positive definite, a condition that 
has been persistently emphasized in the literature. Showing positive definiteness, for a 
radial function, in any d-dimcnsional Euclidean space, is a nontrivial matter. From the 
seminal work of Schoenberg ([5, 25, 26]), we know that p(/3,7) is a completely monotonic 
function if and only if p(f3, 7)(||x|| 2 ) is positive definite on M. d for all d. From the point of 
view of RF theory and modelling, and for the construction of characteristic functions of 




7 



x>0, 



(1.1) 
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spherically symmetric distributed random vectors, it is therefore of interest to determine 
values of the parameters (3, 7 for which p(/3, 7) is completely monotonic. 

Our goal has been to characterize the range of parameters j3, 7 for which p((3,j) is 
completely monotonic and this is done in Theorem 9 below. Necessary conditions for 
complete monotonicity are that /?7 < 1 and j3 < 2, while (3j <1 and j3 < 1 are sufficient. 
When /?7= 1, the functions are not completely monotonic for 1 < (3. However, for (3j < 
1, the picture is more complicated: for 1 < (3 < /3«, we find that p((3, 7) is completely 
monotonic for < 7 < (1 — I (/?))/ ((3 + l(/3)) for a certain continuous and strictly increasing 
function I mapping [1,2] onto [0,2] and such that 1(0*) = 1. 

However, we cannot decide if this is the full range of parameters for which /?(/?, 7) is 
completely monotonic. 

The function I, together with another function c, is used to give a complete character- 
ization of the parameters a,/3 > such that 



and all its positive powers are completely monotonic. We believe that this new class of 
completely monotonic functions may be of some interest in itself. 

The paper is organized in the following way. Section 2 contains basic facts about com- 
pletely monotonic, Sticltjcs, logarithmically completely monotonic and Bernstein func- 
tions in order to make the exposition more readable. In Section 3, we present the auxiliary 
family (1.2) and determine the range of parameters for which f a p is completely mono- 
tonic. In Section 4 we determine the range of parameters for which the functions (1.2) are 
logarithmically completely monotonic and we use these criteria to discuss the parameter 
range for (1.1) to be completely monotonic, hence to obtain permissible correlation func- 
tions associated with a weakly stationary Gaussian RF defined on any Euclidean space. 
In Section 5, we comment on possible comparisons between the Dagum and Cauchy 
models. 

2. Families of completely monotonic functions 

We recall that a function /:]0,oo[— >R is called completely monotonic if it is C°° and 



By the celebrated theorem of Bernstein, the set C of completely monotonic functions 
coincides with the set of Laplace transforms of positive measures fi on [0,oo[ (cf. [28]), 
that is, 



Jo 

where the only requirement on p is that e~ xt is /i-integrable for any x > 0. The set C is 
closed under addition, multiplication and pointwise convergence. 



(_!)«/(«) ( x ) > for x > 0, n = 0, 1, . . . . 



(2.1) 




(2.2) 
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In [10], the authors call a function /: ]0, oo[— >]0, oo[ logarithmically completely mono- 
tonic if it is C°° and 

(-l) n [log/(a;)]^ >0 fori>0 1 n=l,2,... 1 (2.3) 

which are all the conditions of complete monotonicity of log/, except that log / need not 
be non-negative. If we denote the class of logarithmically completely monotonic functions 
by C, then we have / € C if and only if / is a positive C°°-function such that — (log /)' 6 C. 

The functions of class C have been implicitly studied in [1], and Lemma 2.4(h) in that 
paper can be stated as the inclusion £ C C, a fact also established in [9]. 

The class C can be characterized in the following way, established by Horn [15], The- 
orem 4.4. 



Theorem 1. For a function f : ]0, oo[— ► ]0, oof, the following statements are equivalent: 

(i) feC; 

(ii) f a eC for all a>0; 

(iii) V/eC for all n= 1,2,.... 

From Theorem 1, it follows that C is closed under products and pointwise convergence 
since C is. Note that x~ a 6 C for all a > 0. 

We shall now consider an important subclass of C 

Definition 1. A function /:]0, oof— >R is called a Stieltjes function if it is of the form 

/w=0+ r*fm (2 . 4 , 

Jo x + t 

where a > and \x is a positive measure on [0, oof. (The conditions on /i imposed by (2.4) 
can be summarized as J 1/(1 + 1) dfj,(t) <oo.) 

The set S of Stieltjes functions is a convex cone such that S C C. The inclusion is easy 
to verify, but we even have S C £ because for a Stieltjes function /, we have f a £ S C C 
for any < a < 1 (cf. [4] or the survey lecture [5]). 

The formula 

I r°° 

shows that l/x(l + x 2 ) is completely monotonic, but it cannot be a Stieltjes function 
since it has poles at ±i and it follows immediately from (2.4) that a Stieltjes function 
has a holomorphic extension to the cut plane C\] — oo,0]. 

Definition 2. A function /:]0, oo[— ► [0, oof is called a Bernstein function if it is C°° 
and f € C. 
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The set of Bernstein functions is denoted B and is a convex cone closed under pointwise 
convergence. 

Since a Bernstein function is non-negative and increasing, it has a non-negative limit 
/(0+). Integrating the Bernstein representation of the completely monotonic function /' 
gives the integral representation of / € B 

r-OQ 

f(x)=ax + ft + / (l-e- xt )dv(t), (2.6) 
Jo 

where ct,ft> and v, called the Levy measure, is a positive measure on ]0,oo[ satisfying 

OO i 

dv(t) < oo. 

1 + t w 

The following composition result is useful; see [6] or [5]. 
Theorem 2. Let X be any of the sets B,C,C. Then 

feX, geB^fogeX. 



3. The auxiliary family F 

The main result of this paper, regarding the Dagum family as described in equation 
(1.1), is obtained by making use of the following family of functions, indexed by the two 
parameters a, ft > 0: 

F:={f aAx) = ^^ yX >0}. (3.1) 

We want to discuss the range of parameters for which these functions are completely 
monotonic. 

Since a completely monotonic function has, by (2.2), a holomorphic extension to the 
half-plane SRz > and since z 13 = — 1 in this half-plane if (3 > 2, it follows that (3 < 2 is a 
necessary condition for f a ,p G C. 

Theorem 3. The following statements hold: 

(i) f a , p eC for a> 0,0 <0<1; 

(ii) f a ,peC for a>/3/2,l</3<2. 

Proof. Clearly, x~ a € C for a > 0. Since x 13 € B for < /3 < 1 and the reciprocal of a non- 
zero Bernstein function is completely monotonic (cf. [6], page 66) we get 1/(1 + x 13 ) € C 
for these values of ft. Using that C is closed under products, we see that (i) holds. 

Composing the completely monotonic function (2.5) with the Bernstein function 
a^O < ft < 1, we get fp,2/3 G C by Theorem 2. Multiplying f/3,20 by x~ a , we get (ii). □ 
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We define the function c : [1, 2] — * [0, oo[ as follows: 

c(/3)=inf{a>0|/ Q , /3 eC}. (3.2) 

By Theorem 3(h), we get c(/3) < /3/2 < 1 and since C is closed under pointwise conver- 
gence, we get f c (p),0 G C, hence 

f a , P eC^a>c(p), l</3<2. (3.3) 

While we know from Theorem 3 that c(l) = and c(2) < 1, we next prove that c(2) = 1. 

Lemma 1. 

/ a ,2 6C^a>l. 

Proof. We have to show that f a ,2 C for < a < 1. This is clear for a = since 

poo 



1 



1 + x 2 

For < a, we define the function n a by 



e sinidi. 



o 



*a(t) = S-r, *>0, (3.4) 
r(a) 



and it is easy to see that 

1 



e -xt t a-l dt 



r(a) y 

that is, that x~ a is the Laplace transform of K a , hence 

i- 1 



(3.5) 



f a ,2( x )= e- xt r) a (t)dt, r )a (t) = —^ (t-s) a - 1 sinsds. 
Jo T ( a ) Jo 

The assertion follows if we prove that r] a (2n) < for < a < 1 . We find 

/•27T |>7T 

r(a)77a(27t) = / (27t-s) a - 1 sinsds= / ((2tt- s)"" 1 - (tt - s)"" 1 ) sinsds < 
Jo Jo 

because (2n— s) Q_1 <(n — s)" -1 and sins > for < s < 71. □ 

Formula (2.5) can be generalized. For 1 < (3 < 2, we want to find the function <pp such 
that for x > 0, 

oo 

e- xt ipf3{t)dt. (3.6) 



1 + xP 
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Theorem 4. Let 1 < (3 < 2. The function ipp of (3.6) is given by 
sin(/37r) f°° e- ts s ds 



o 1 + 2sP cos(/?7t) + s 2 ? 







71 



( 71 



— exp t cos [ —— ) i cos — + 1 sin — 



\0 



t>0. 



(3.7) 



Proof. Fix (el. The function f(z) = c t2 /(l + z@) is meromorphic in the cut plane 
A = C\] — oo,0] with simple poles at z = exp(±i(7t//3)). Defining 



00 pity 



1 f°° 



e lty dy 
7^ 



(3.8) 



wc get (pp (t) = for t < by integrating / along the half-circle in the right half-plane 
with diameter [—iR, iR] and letting R — > oo. For t > 0, we integrate / along the contour 
consisting of 



-LR,ii?], z = Re w , 0e 
z = re ±i(n-e) 



71 

2^ 



U 



7T 

-n + s,-- 



■e[0,R], 



where R > 1,0 < e < tc(1 — 1/0), and next we let R — > oo and e — > 0. Since the residues 
of / at the poles are 

Res(/,exp(±i(7t//?))) = -i c ^/P) exp^e^ 7 ^), 

we get the expression (3.7) for <fp{t). By the Fourier inversion theorem, 1/(1 + (iy) ) is 
the Fourier transform of ipp defined by (3.8) and this shows (3.6). □ 

Remark 1. It is clear that (3.6) holds for (3 — 1 with ipi(t) = cxp(— t) and for (3 = 2 with 
f2(t) = sint. The latter also follows immediately from (3.7) by letting (3^2. However, 
it is not directly possible to let (3 — > 1 in this formula. After an integration by parts using 



/3s - 1 ds 



1 



1 + 2 cos(/37t)s/ 3 + s 2 f j sm((3n) 



arctan 



s> 3 + cos((3n) 
sm((3n) 



wc find 



(3n 



arctan ■ 



■ cos((3n) 



sin(/?7t) 



71 



(1 -~ts)e~ ts ds 



— exp t cos -r cos — + t sin 



i>0, 



(3.9) 



and splitting the integral into integrals from to 1 and from 1 to oo, it is possible to let 
/3 — ► 1 and to obtain exp(— t) in the limit. 
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Generalizing (3.5), we have 

f°° 1 /"* 

faA x )= e-^ffc^Jdt, W (i) = -— / (t- s ) Q -V/3(s)d s (3.10) 
Jo r (a) Jo 

and f a . p £ C if and only if rj a ^(t) > for t £ [0, oo[. 

Remark 2. We have c(/3) > for 1 < (3 < 2. In fact, since ifp(t) < in the t-intervals 
where cos(S+tsinS) > 0, we see by continuity that if wc fix such a t, then also r) at p(t) < 
for a sufficiently close to 0. Here, we use the fact that lim Q ^ K a = <$o vaguely (cf. [6], 
Proposition 14.26). 

Theorem 5. The function c: [1,2] — ► [0, 1] is strictly increasing, lower semi-continuous 
and continuous from the left. 

Proof. Fix 1 < Pi < 02 < 2. We compose f c (0 2 ),0 2 ^ C with the function x^ 1 ^' 2 and get 

by Theorem 2, hence c(/?i) < /3ic(/3 2 ))//?2 < c(/3 2 ) as asserted. Fix 1 < [3 < 2 and < a < 
c(/?o)- Then / a!l g ^ C and using the fact that C is closed under pointwise convergence on 
]0, oof, wc see that f a ,p ^ C for (3 in a sufficiently small neighborhood V of (3q relative to 
]1, 2], hence c((3) > a for /3 £ V, but this proves that c is lower semi-continuous at /3 - An 
increasing and lower semi-continuous function is continuous from the left. □ 

Remark 3. We conjecture that c is continuous but have not been able to verify this. 



4. The Dagum family and the class C 

In this section, we inspect the range of parameters for which the Dagum class, defined 
in equation (1.1), is completely monotonic on the positive real line. 

Let us begin by showing the range of parameters a, (3 for which f a jj £ C. 

Theorem 6. The following statements hold: 

(i) f a , ££ for a >0,0<P<1; 

(ii) / Q , 2 e£^a>2. 

Proof, (i) By Theorem l(iii), it suffices to prove that (x a (l + a;' 3 )) -1 /" £ C for n = 
1,2,..., but this holds because (1 + a;' 9 ) 1 /" £ B by Theor em 2 and the reciprocal of a 
Bernstein function is completely monotonic. 
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(ii) We shall determine the parameters a > such that — (log/c^)' € C, that is, such 
that 

6C. 



:r 1 + x 2 

However, this function is indeed the Laplace transform of a + 2 cost, which is non-negative 
on [0, oo [ if and only if a > 2. □ 

We now define the function I : [1, 2] — > [0, oo[ by 

lC9)=inf{a>0|/ a ,/jG£}. (4.1) 

Using the fact that C is closed under pointwise convergence, we get fi(p)^ G £■ and 
observing that if f at p £ C, then / Q ',/3 € £ for any a' > a, we get that 

f a ,p££^a>l([3), l</3<2. (4.2) 

From Theorem 6, we know that 1(1) = 0, 2(2) = 2 and, clearly, c(f3) < l((3). 
For 1 < j3 < 2, we introduce the function 

l/»/j(t)= / ^(s)ds, i>0. (4.3) 
Jo 

Since it is equal to n\ * tpp, we know that 

1 



x(l+xP) 



= UiPp](x), (4.4) 



which is completely monotonic because c(f3) < 1, hence i^p(t) > for i > 0. From formula 
(3.7), we find (l</3<2) 

Mt)=Pf3(t)+Tp(t), (4.5) 

with 

pp(t) = 1 - -^expficos^ cos sin ^ (4.6) 

and 

sin^Tt) r e^V^dg 
™ 7t io l + 2^cos(/37t) + S 2 / 3 1 j 

because ^ = ipp and ^(O) = 0. Note that rg(i) e C with ^(0) = 2//3 — 1 and 

W (dn^) =1 + | CXP ( 7TCOt ?) >1 ' 
hence tpp(t) nas a global maximum which is greater than 1 on [0, oof. 
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Substituting s 13 = —usm((3n) — cos(/37t) in (4.7), we get 

= i r» eM -t(-uMM-cos(^) ^ 

P /W-cot(^) l + « 2 

Let U = UCB([0,oo[) denote the Banach space of uniformly continuous and bounded 
real-valued functions on [0, oo[ with the uniform norm. For each 1 < (3 < 2, it is clear that 
ipp € U since 

/>oo 

lim ^(i) = / (pp(s)ds = l. 



Theorem 7. The function (3 — > ipp from ]1, 2[ to U is continuous. Furthermore, 



and 



limV^(t)=V , i(*) = l-e _ * mt7 (4.9) 



lim ipp{t) = tp2(t) = 1 — cost locally uniformly for t G [0, oof. (4-10) 

/3 >2 



Proof. The function pp(t) is well defined and continuous for ((3,t) £ [1,2] x [0,oo[ with 

pi{t) = 1 - 2c _t , p 2 (t) = 1 - cost. 

It is easy to see that (3 —> pp is continuous on [1,2[ with values in the Banach space 
U, while pp(t) — ► P2(t) locally uniformly on [0,oo[ when (3 — >2. (We cannot have global 
uniform convergence since /9g tends to 1 at infinity while p^ is oscillating.) By (4.8), it is 
easy to see that (3 — * Tp is continuous from ]1, 2[ to U and that 

lim ra(t) = e _t , limT/3(i) = in [/, 

/9— »1 /3^2 

that is, uniformly for £ G [0,oo[. □ 
Corollary 1. TTie function f3 —> ^(/3) — max t >o ipp (t) is continuous on [1,2]. 
Proof. We know that *(1) = 1 and that *(/?)> 1 for 1< [3 < 2. Furthermore, 

*(/?) < nm^(i) + maxr^) <l + | + |- l = l 

If F : X — > [/ is a continuous mapping of a topological space X with values in the Banach 
space U, then supF:X^M is also continuous. Therefore, \& is continuous on [1,2[. By 
(4.10), we get 

max il>d(t) — > max "02 (t) = 2. 

0<t<27t ' 0<t<27t 



1144 



C. Berg, J. Mateu and E. Porcu 



hence liminf ( g^2 ^{P) > 2. On the other hand, limsup^_, 2 ^W) ^ lim J g_ > 2 4 = 2 and a 
combination of these results gives lim^ 2 ^{(3) = 2. □ 

Theorem 8. T/ie function I is given by 

l(0)=0(*(J3)-l), l</3<2, (4.11) 
and it is continuous and strictly increasing. Furthermore, ^ is increasing. 

Proof. We have 

-(log /„,„)'(*) = — - n f = L[(a +/3) - /%](*) 

and therefore / QiJ g S £ if and only if ^(i) <l + a/(3 for £ > 0, which holds if and only if 

a >/?(*(£) -1), 

giving the formula (4.11). The continuity of Z follows from Corollary 1. 

We prove that l((3)/(3 is increasing, which implies that Z(/3) is strictly increasing and 
that <I> is increasing. Let 1 < /3i < fa < 2. We compose fic/3 2 ),p 2 € £ with the function 

a;/3i /ft and get 

by Theorem 2, hence l(fa) < fal{fa)/fa < l(fa) as asserted. □ 

Using that the function Z is continuous and strictly increasing such that Z(l) = 0, 1(2) = 
2, the following definition makes sense. 

Definition 3. Let fa denote the unique number in ]1,2[ such that l(fa) = 1. 

The number fa is where the graphs of ^{(3) and 1 + 1/(3 intersect, and from Figure 
1 (produced using Mathcmatica), we estimate that fa is approximately 1.74. For fa < 
(3 < 2, we have Z(/3) > c(/3). It is reasonable to believe that this inequality holds also for 
1 < (3 < fa. In any case, for 1 < (3 < 2, we have 

f a , eC\C for c(/3) < a < Z(/3) (4.12) 

and 

for 7 >0,«>Z(/3). (4.13) 

We are now able to state the main result of this paper, regarding the Dagum family 
p(j3, 7) as defined in (1.1). 
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Theorem 9. Let /3,7> 0. We then have 

T /37-l 

P (P, 7 )eC^ {1 + x0)y+1 eC. (4.14) 

// these conditions hold, then (3^ < 1 and /3 < 2. Ira i/ie opposite direction, we have 

(i) if Pi < I and f3<l then p(/3, 7) e C ; 

(ii) «/ /37 = 1, then 

P (f3,l/(3)€C&/3<1; (4.15) 

(iii) z/ /?7 < 1 and 1 </3</3*, i/ien p{(3,~f) € C i/, in addition, 

Proof. Using the fact that p(/3, 7) is a positive C°°-function, it is completely monotonic 
if and only if — p(/3, 7)' <G C, and a small calculation leads to the right-hand side of (4.14). 

If this function is completely monotonic, the value at x — has to be strictly positive 
or infinite and hence /?7 < 1. Furthermore, (3 < 2, otherwise the function has poles in the 
right half-plane and then it cannot be completely monotonic. 

If /?7<1 and (3<1, then 

eC (4.17) 



(1 + 0^)7+1 Vz (1-/37)/(1+7) (l + ^) 




~i 1 1 1 r 

1.2 1.4 1.6 i.b 2.0 
I ;■■!.-! 




-1 1 1 1 r 

1.2 1.4 1.6 1.8 20 
I ■■■ ■ '<■■< 




1.2 14 1.6 1.8 2.0 



Figure 1. From upper left to lower right: max(\E>(/3)) vs. /3; 1 + 1//3 vs. (3; intersection of the 
first two graphs; l(/3) vs. j3. 
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because f a ^ £ C for any a > 0, which shows (i). 

To see (ii), it suffices, for c > 0, to prove the following implication: 

f§ t p€C=^0<l. (4.18) 

Assuming /g ^ G C, we know that (/q »)" > 0, hence 

PdP-*{\ + x' 3 )- c - 2 [2; & (/3c + 1) + 1 -0\ > 0, x > 0, 



and therefore 



lim [a^(/3c+l) + l-/?]>0. 

x— >0+ 



This shows that /3 < 1 . 

The assertion (iii) also follows from formula (4.17) because 1(0) < 1 for 1 < (3 < /3* and 

Theorem 9 shows the range of completely monotonic permissibility for the Dagum 
class. Under these conditions, we have that this function is the correlation function of an 
isotropic and Gaussian RF. Rephrasing this, the result shows that there exists a (wide 
sense) stationary Gaussian RF defined on any d-dimensional Euclidean space, whose cor- 
relation function is defined by the Dagum class. The same result admits an interpretation 
in terms of characteristic functions of positive spherically symmetric distributed random 
vectors of M. d for all d £ Z + . 

Remark 4- For j3 = 2, results about the complete monotonicity of 

1 



x a (l + x 2 ) x 

have been obtained in a number of papers. We mention the following results: 

(i) g a ,x G C for a > A > 1 (cf. [11]); 

(ii) g a , x G C for a > 2A > (cf. [2, 3]); 

(iii) . 9 i,A G C for < A < 1 (cf. [19]); 

(iv) <? Q: a ^ C for < a < A (cf. [19, 29]) - the simple proof in [29], Lemma 8, carries 
over to this situation; 

(v) 5q , q i C for < a < 1 (cf. [19]). 



5. Comparisons between models and some conclusions 

In this paper, we have presented some results for the range of parameters allowing the 
complete monotonicity of the Dagum function. Some comments arc in order. 
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A first consideration is that there is not, to the knowledge of the authors, any compar- 
ative study available in the literature regarding the comparison of the Cauchy and the 
Matern models. The reason for this is that the Matern class possesses an explicit analyt- 
ical form for the associated spectral density and this allows some important properties of 
the associated random field (RF) to be known. Such properties cannot be established in 
the case of the Cauchy model, as the Fourier pair can be calculated only for some partic- 
ular parameter setting. The same holds for the Dagum class. In addition, the Matern and 
Cauchy correlation functions represent the most widely used families for geostatistical 
modelling. 

A natural comparison arises between the Dagum and Cauchy correlations, being both 
models decouplers of the fractal dimension and Hurst effect. With this purpose, in a 
recent paper ([17]) an extensive simulation study is performed in order to highlight the 
attitude of the Dagum correlation function with respect to decoupling. In this case, the 
authors use the convenient paramctrization 



One can verify that sufficient conditions for the permissibility of this function on any 
d-dimensional Euclidean space are 7 € (0, 2] and e < 7. Recall that the Cauchy model 
has expression 



Observe that the parameters identifying the fractal dimensions associated, respectively, 
with Dagum and Cauchy classes, are e, 8 <G (0, 1], while those identifying the Hurst effect 
are (respectively) S, r\ <E (0, 1). From the results contained therein, it is quite evident that, 
when simulating profiles (i.e., in the one-dimensional case) the Cauchy class shows better 
performances in comparison with the Dagum class in terms of decoupling. 

The dilemma arises in two- or higher-dimensional spaces, where it is far from clear 
which of the two models performs better. Reaching conclusions, in this case, is not a 
trivial matter, as several variables can influence the performance of two decouplers in 
two- or higher-dimensional spaces and the nature of these variables is very often of a 
non-statistical nature. This is verified, for instance, in [20, 21]. Thus, everything seems 
to converge to some open problems whose solution could be found through heuristic 
arguments and deeper simulation studies, taking into account the influence of some real 
variables on the behavior of the decouplers in bidimensional spaces. 

In addition, and as far as estimation of the parameters is concerned, we performed 
several simulation studies using likelihood approximations such as the weighted composite 
likelihood, proposed by [7], and for both models the results are very promising in the 
sense that, in both cases, the parameters identifying fractal dimension and Hurst effect 
are estimated with a high level of accuracy. 

Thus, the Dagum model behaves similarly to the Cauchy model in terms of decoupling 
and estimation properties. In addition, the Dagum function is actually opening several 
avenues of research in both applied and theoretical studies. For the latter, initial evidence 





P (t) = (i+t e r^ e , 



t:= j|h|| 6 1^ Be (0,2],jj>0. 
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of this fact can be found in [14] , where the Cauchy and Dagum functions are studied in 
terms of extremal rays of positive definite functions. As an example of the former, the 
Dagum model is lately being considered in material science and mechanical engineering 
(see, e.g., [21, 22]), where practical analysis in this field is performed using the Dagum 
function. These considerations reinforce the added value of the Dagum model, which can 
be considered a good alternative to such other important published models. 
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